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ABSTRACT 

We report observations of the nova RS Ophiuchi (RS Oph) using the Keck 
Interferometer Nuller (KIN), approximately 3.8 days following the most recent 
outburst that occurred on 2006 February 12. These observations represent the 
first scientific results from the KIN, which operates in N-band from 8 to 12.5 
/im in a nulling mode. The nulling technique is the sparse aperture equivalent of 
the conventional coronagraphic technique used in filled aperture telescopes. In 
this mode the stellar light itself is suppressed by a destructive fringe, effectively 
enhancing the contrast of the circumstellar material located near the star. By 
fitting the unique KIN data, we have obtained an angular size of the mid-infrared 
continuum of 6.2, 4.0, or 5.4 mas for a disk profile, gaussian profile (FWHM), 
and shell profile respectively. The data show evidence of enhanced neutral atomic 
hydrogen emission and atomic metals including silicon located in the inner spatial 
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regime near the white dwarf (WD) relative to the outer regime. There are also 
nebular emission lines and evidence of hot silicate dust in the outer spatial region, 
centered at ~ 17 AU from the WD, that are not found in the inner regime. 
Our evidence suggests that these features have been excited by the nova flash 
in the outer spatial regime before the blast wave reached these regions. These 
identifications support a model in which the dust appears to be present between 
outbursts and is not created during the outburst event. We further discuss the 
present results in terms of a unifying model of the system that includes an increase 
in density in the plane of the orbit of the two stars created by a spiral shock wave 
caused by the motion of the stars through the cool wind of the red giant star. 
These data show the power and potential of the nulling technique which has 
been developed for the detection of Earth-like planets around nearby stars for 
the Terrestrial Planet Finder Mission and Darwin missions. 

Subject headings: novae: cataclysmic variables — stars: dwarf novae — stars: 
individual (RS Ophiuchi, HD 1622 14) — techniques: interferometric — techniques: 
high angular resolution 



1. INTRODUCTION 

Classical novae (CN) are categorized as cataclysmic variable stars that have had only 
one observed outburst - an occ urrence typified by a Johnson V-band brightening of between 



six and nineteen magnitudes (jWarner I Il995l ). These eruptions are well-modeled as ther- 



monuclear runaways (TNR) of hydrogen-rich material on the surface of white dwarf (WD) 
primary stars that, importantly, remain intact after the event. Current theory tells us that 
CN can be modeled as binary systems in which a lower-mass companion - the secondary - 
orbits the WD primary such that the rate of mass transfer giving rise to the observed erup- 
tion is very low. Recurrent novae (RN) are a related class of CN which have been observed 
to have more than one eruption. Like CN, RN events are well-represented as surface TNR 
on WD primary stars in a binary system, but are thought to have much higher mass transfer 
rates commensurate with their greater eruption frequency. There are two types of systems 
that produce recurrent novae - CVs, in which the WD accretes from a main sequence star 
that orbits the WD on a time scale of hours, and symbiotic stars, in which the WD accretes 
from a red-giant companion that orbits that WD on a time scale of years. 

CN and RN produce a few specific elemental isotopes by the entrainment of metal- 
enriched surface layers of the WD primary during unbound TNR outer-shell fusion reactions. 
In contrast, type la supernovae produce most of the elements heavier than helium in the 
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Universe through fusion reactions leading to the complete destruction of their WD primary. 
Some theoretical models indicate that RN could be a type of progenitor system for super- 
novae. Importantly, these theories are predicated on two critical factors: 1.) the system 
primary must be a compact carbon-oxygen core supported solely by electron degeneracy 
pressure and 2.) there must be some mechanism to allow the WD mass to increase secularly 
towards the Chandrasekhar limit. 

The nova RS Q ph has undergone six recor ded epi sodic outbursts o f irregu lar interval in 

1898 jFleminsJE^oj^ pg jAdams fc JovlEoSsh . 1958 (|Wallersteinlll958h . 1967 JSarbon. Mammano fc Rosi 

1969h. 1985 JMorrisonl Egsi) and now 2006. There are al so two possible outbursts in 1907 

( jSchaeffer 2004) and 1945 ( Oppenheimer &: Mattel 1993). All outbursts have shown very 

similar light curves. This system is a si ngle-line binary, symbiotic with a red giant sec- 

ondary characterized as K5.7 ± 0.4 I-II (IKenyon fc Fernandez- Castrol 119871 ) to a K7 III 

(IMiirset fc Schmidlll999l ) in quiescence and a white dwarf primary in a 455.72±0.83 day orbit 

about their comrn o n center of mass as meas ured using single-line radial velocity techniques 

(IFekel et al.ll2000l ). lOppenheimer fc Mattel! examined all the outbursts and found that V 

band luminosity of RS Oph decreased averaged 0.09 magnitudes/ day for the first 43 days 

after outburst. A 2- magnitude drop would then require on average 22 days, establ ishing RS 

Oph as a fast nova based on the classification system of iPayne-GaposchkinI (119571 ). 



The most recent outburst of the nova RS Oph was discovered at an esti mated V-ban d 
magnitude of 4.5 by H. Narumi of Ehime, Japan on 2006 February 12.829 UT (jNarumill2006l ). 
This is 0.4 mag brighter than its historical average AAVSO V-band peak magnitude so it is 
reasonable to take Feb 12.829 (JD 2453779.329) as day zero. The speed of an outburst is 
characterized by its t2 and ts times which are the intervals in days from the visible maximum 
until the system has dimmed by 2 and 3 magnitudes, respectively. For this outburst the ^2 
and ^3 times are 4.8 and 10.2 days, respectively. 

The distance to the RS Oph system is of impo rtance to the interfer ometry community 
as it effects interpretation of astrometric data (cf. iMonnier et al.l (120061 )). There has been 



a good dea l of disagreement in th e literature with a s urprisingly broa d range, frona as near 
as 0.4 kpc JHachisu fc KatolboOlh to as far as 5.8 kpc jPottaSIlSS). ISarrv et al.l J2008ah 



have recently undertaken a thorough review of the various techniques that have been used 
to derive a distance to RS Oph and obtain a distance of 1.41q'2 kpc. It is this value that we 
adopt for astrometric calculations in this paper. 



The structure of this paper is as follows. We report high-resolution N-band observations 
of RS Oph using the nulling mode of the 85 meter baseline Keck Interferometer, beginning 
with a discussion of the nulling mode itself in Section 2. We discuss the observations in 
Section 3, and the data and analysis in Section 4. In Section 5 we introduce a new physical 
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model of the system, which unifies many of the observations into a coherent framework. The 
results of this paper and those of other recent observations of RS Oph are discussed in the 
context of this model in Section 6. Finally, Section 7 contains a summary of our major 
results and conclusions. 



2. THE KECK INTERFEROMETER NULLER 

The KIN is designed to detect faint emission due, e.g., to an optically-thin dust envelope, 
at small angular distances from a bright central star (Serabyn, Colavita and Beichman 2000). 
Its operation differs from a more common fringe scanning optical interferometer in that a 
nulling stage precedes the fringe scanning stage (Serabyn et al. 2004, 2005, 2006; Colavita et 
al. 2006). The basic measurement thus remains the fringe amplitude, but both the meaning 
of the fringe signal in relation to the source and the processing of the fringe information 
differ from the normal case of a standard visibility measurement. Here we provide only a 
brief description of the measurement process, because this has been and will be described in 
depth elsewhere (Serabyn et al. 2004, 2005, 2006, Colavita et al. 2006; Serabyn et al. 2008, 
in prep). 

To remove both the stellar signal and the thermal background in the MIR, a two stage 
interferometer has been developed. To implement this approach, each Keck telescope is 
first split into two half-apertures, to generate a total of four collecting sub-apertures. The 
starlight is then first nulled on the two long (85 m) parallel baselines between corresponding 
Keck subapertures. This generates the familiar sinusoidal fringe pattern on the sky (Fig. [1]), 
except that the central dark fringe on the star is achromatic, and fixed on the star, to achieve 
deep and stable rejection of the starlight. After the nulling stage, the residual, non-nulled 
light making it through the first stage fringe pattern is measured by a fringe scan in a second 
stage combiner, the "cross-combiner", which combines the light across the Keck apertures 
(~ 4 m baseline). Thus, what is measured is the fringe amplitude of this "nulled source 
brightness distribution" (Serabyn et al. 2008, in prep.) This quantity is then normalized 
by the total signal. This is measured by moving the nullers to the constructive phase, and 
again scanning the cross-combiners. The basic measured quantity, the null depth, A^, is then 
the ratio of the signals with the star in the destructive and constructive states. N is related 
to the classical interferometer visibility V = {Imax — I mm) / {I max + -^mm), the modulus of the 
complex visibility V , by a simple formula: 



(1) 
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The long nulling baselines produce fringes spaced at about 23.5 mas at 10 /im, while 
the short baseline produces fringes spaced at ~ 400 mas, which is similar to the size of the 
primary beam and is assumed to be large compared to the extent of the target object. 
Modulating its phase therefore modulates the transmitted flux of the entire astronomical 
source, as modified by the sinusoidal nuller fringes, and so the amplitude of this modulated 
signal gives the flux that passes through the fringe transmission pattern produced by the 
long baseline nuUers. 



3. OBSERVATIONS 

We observed the nova about an hour angle of about -2.0 on the Keck Interferometer in 
nulling mode on 2006 Feb 16 with a total of three observations between day 3.831 and 3.846 
post-outburst bracketed with observations of two calibrators stars, p Boo and x UMa. Data 
were obtained at N-band (8 to 12.5 /im) through both ports of the KALI spectrograph with 
the gating of data for long baseline phase delay and group delay turned on. The Infrared 
Astronomical Satellite (IRAS) Low Resolution Spectrometer (LRS) spectra for the calibrator 
stars were fiux-scaled according to the broadband IRAS 12 yum fluxes. The nova data were 
flux calibrated and telluric features were removed using calibrator data. These calibrators 
are well matched to the target flux and their size and point- symmetry are well known. A 
journal of observations is presented in Table [TJ 

Our data analysis involves removing biases and coherently demodulating the short- 
baseline fringe with the long-baseline fringe tuned to alternate between constructive and 
destructive phases, combining the results of many measurements to improve the sensitivity, 
and estimating the part of the null leakage signal that is associated with the finite angular 
size of the central star. Comparison of the results of null measurements on science target 
and calibrator stars permits the instrumental leakage - the "system null leakage" - to be 
removed and the off-axis light to be measured. 

Sourc es of noise in the me asurements made by this instrument have been well described 
elsewhere ( iKoresko et al. I l2006h . however, we outline them here for reference. The null leakage 



and intensity spectra include contributions from the astrophysical size of the object, phase 
and amplitude imbalances, wavefront error, beamtrain vibration, pupil polarization rotation, 
and pupil overlap mismatch. There are also biases, mostly eliminated by use of sky frames, in 
the calculated fringe quadratures, caused mainly by thermal background modulation due to 
residual movement of the mirrors used to shutter the combiner inputs for the long baselines. 
Another source of error is the KALI spectrometer channel bandpass, which is large enough 
to produce a significant mismatch at some wavelengths between the center wavelength and 
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the short-baseline stroke OPD. This effect, termed warping, distorts the quadratures and 
is corrected by a mathematical dewarping step accomplished during calibration. There is 
also the effect of the partial resolution of extended structure on the long baselines which 
will cause the flux to be undercounted by some amount depending on the spatial extent 
and distribution of the emitting region. Compared to the error bars this is a rather small 
effect for most normal stars, and is unlikely to have much influence on the actual spectrum, 
though, unless the emission lines are coming from a very extended shell approaching 25 mas 
in angular size. We do not expect this to be the case for nova RS Oph at day 3.8. 

The most important contributor to measurement noise is sky and instrument drifts be- 
tween target and calibrator. Other potential sources of noise include the difference between 
the band center wavelength for the interferometer and that of the IRAS LRS calibrator, 
undercounting of stellar flux resulting from glitches in the short-baseline phase tracking, and 
the chromaticity of the first maximum of the long-baseline fringe. None of these are signifi- 
cant for the following reasons. First, the calibration is based not on broadband photometry 
but on the IRAS LRS spectrum. Second, short-baseline tracking glitches happen nearly as 
frequently on target observations as on calibrator observations so they should have minimal 
effect on overall calibrated flux. Third, the effect of chromaticity should be negligible because 
fringe detection is done on a per-spectral-channel basis. As a result, it is only affected by 
dispersion within the individual KALI spectral channels, which are about 0.3 fj^m in width. 



4. DATA & ANALYSIS 

We developed a mathematical solution and software suite to model the observatory and 
source brightness distribution. We used this suite to conduct an exhaustive grid search and 
to generate a monte carlo confidence interval analysis of solution spaces of these models. We 
explored three types of models for the source surface brightness distribution; Gaussian, disk, 
and shell. Limited (u,v) coverage permitted only rotationally-symmetric models with two 
parameters - size and flux. We used minimization to obtain the best fit models for both 
the inner and outer spatial regimes simultaneously. Table [2] displays size measurements, flux 
values, and one-o" confidence interval values. The error bars have been increased slightly 
beyond the values to incorporate the effect of an adopted 0.005 one-sigma systematic 
error which would be correlated among measurements at different wavelengths. 

The measurements made with the two KALI ports are somewhat independent - the data 
they produce are combined for purposes of fringe tracking, but not for data reduction. The 
system null and the final calibrated leakage are computed separately for the two ports. The 
apparent inconsistency detected between the ports is the result of optical alignment drifts 
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at the time of the measurement. In particular, the last calibrator measurement showed a 
sudden change in the system null for Port 1, while for all the other calibrator measurements 
the system nulls were stable. We therefore compared our source brightness distribution 
models against KALI port 2 data alone. For the best-fit models in Table [2] we used Spitzer 
spectra to identify and remove emission features centered at 8.7, 9.4, 10.4, 11.4, and 12.5 /zm 
in the KIN inner and 8.9, 9.8, and 11.4 /im in the KIN outer spatial regime. We removed 
the emission feature data because our intention was to model the continuum. 

Figure [2] shows two sets of 8 — 12.5 /im spectra of RS Oph on day 3.8 post-outburst. The 
upper plot shows the outer spatial regime, which is the dimensionless null leakage spectrum, 
i.e., the intensity of light remaining after destructive interference divided by the intensity 
spectrum, plotted against wavelength in microns. The lower plot is the intensity spectrum, 
which is light principally from the inner 25 mas centered on the source brightness distribution 
orthogonal to the Keck Interferometer baseline direction - 38 degrees East of North. The null 
leakage spectrum may be broadly described as a distribution that drops monotonically with 
increasing wavelength overlaid with wide, emission-like features. The intensity spectrum, 
with an average flux density of about 22 Jy over the instrument spectral range, may be 
similarly described but with a continuum that has a saddle shape with a distinct rise at 
each end. The underlying shape curves upward for wavelengths shorter than approximately 
9.7 nm and longer than 12.1 /xm. Overlaid on each of these are traces representing simple 
models of source brightness distributions fit to the data, described above. 

Figure [3] shows the in ner and outer spat ial regimes from the KIN data together with 



Spitzer data from day 63 (lEvans et al.l l2007bl ) . The absolute fluxes of these data are un- 



sealed and given with a broken ordinate axis for clarity. The outer KIN spectrum flux has 
been multiplied by 2 to correct for the transmission through the fringe pattern for extended 
sources. In efforts to identify the sources of these emission features we took the high reso- 
lution Spitzer spectra and, using boxcar averaging, re-binned the data until it and the KIN 
data had equivalent resolution. The binned Spitzer spectrum is quite similar in character to 
the KIN spectra. Importantly, the sum of the inner and outer spatial regime KIN spectra is 
nearly identical to the binned Spitzer spectrum with the exception that the absolute scale 
magnitude is, on day 3.8, about an order of magnitude greater than that of the Spitzer spec- 
trum. This is as expected because the measured flux drops with time after the peak due to 
cooling and the summed inner and outer spatial regime data should very nearly reproduce 
the transmission of a filled aperture telescope of equivalent diameter. 

The wavelength range sampled by the KALI spectrometer, 8 - 12.5 fim, covers many 
important discrete transitions including molecular rotation- vibration, atomic fine structure, 
and electronic transitions of atoms, molecules, and ions. This range also samples several 
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important transitions in solids such as silicates found in dust and polycyclic aromatic hy- 
drocarbons (PAHs). With the exception noted below, spectral features in the KIN data 
are clearly not resolved by the instrument and are often, in the associated Spitzer spectrum, 
Doppler broadened and blended. Also, because the Spitzer spectra, were not taken contempo- 
raneously with the KIN spectra and because of the transient nature of the RN, identification 
of KIN features with Spitzer lines is necessarily tentative. 

Figure H] shows Spitzer spectra on days 63, 73 and 209 after peak V-band brightness. 
The continuum decreases monotonically with time. The continuum values are approximately 
1.4, 1.1, and 0.25 Jy for April 14, 26 and September 9, respectively. The first two spectra on 
April 14 and 16 show strong atomic lines with no obvious evidence of dust emission. However 
it is expected that thermal bremsstrahlung from the central source will overwhelm any other 
faint sources. Additionally, there are several narrow emission features in this latter spectrum 
which appear to be similar to those in the earlier spectra. In contrast, the spectrum taken 
on September 9 shows a distinct broad emisison feature from 8.9 to 14.3 /im peaking at 10.1 
fim. We proceed on the assumption that this is emission from dust in the vicinity of the 
nova. 

Referring again to the KIN spectrum in Figure [3l strong continuum radiation is apparent 
in both inner and outer spatial regimes. In contrast, while the contin uum was still clearl y 
visible in J, H, I, K bands on February 24 and detectable on April 9 (lEvans et al.ll2007al ). 
it has subsided by day 63 in Figure Bl and F igure HI Comparing the spectra from RS Oph 
to those from V1187 Sco (ILynch et al.ll2006l ) note that the continuum given for V1187 Sco 
is non-Planckian showing an excess longwards of 9 /im and is strongly red as compared to 
a F5V Kurucz spectrum. The V1187 Sco and RS Oph continua have slopes that agree to 
within 10%. While both free-bound and free-free transition processes lead to emission of 
continuum radiation, in the MIR spectral range thermal bremsstrahlung free-free emission 
dominates. We attribute the drop in continuum radiation to the transition to line-emission 
cooling mechanisms. Also, by the time Spitzer data were taken, the object had become less 
dense and so the emission coefficient for thermal bremsstrahlung (proportional to number 
density of protons an d electrons) had dr opped considerably. The continuum emission is 
described in detail by iBarry et al.l (l2008bl ). 



Table [3] gives identification of all narrow features in Figure H] where it is possible to do 
so, as follows. We generated line lists of atomic species assuming that all ionization stages 
of all elemen ts were possib le for the three Spitzer spectra. We assumed standard cosmic 
abundances (jGrevessdll984l ). After continuum normalization, each emission line was fitted 
with a Gaussian and, where necessary, de-blended using IRAF. Our fitted emission features 
were compared with the Spitzer list. 
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Table H] gives our identification of particular emission sources with the continuum- 
normalized spectral features in the KIN inner and outer spectra. The center wavelength 
of each of the broad features is listed in the first column, while the second column displays a 
width for each feature defined as the cutoff wavelength at the intersection of the the feature 
and the unity continuum line - the full-width zero intensity (FWZI) level of the feature. The 
fiux and the one-o" uncertainty contained in that feature through measurement of the total 
area under it and above the unity continuum level is in the third column. The fourth column 
shows the particular atomic species found in the Spitzer spectrum with each KIN spectral 
feature. Because Spitzer lacks the spatial resolution to discriminate the inner from the outer 
regions of the nova and the KIN lacks the spectral resolution to discriminate among atomic 
species we assumed that features in the KIN spectra would most reasonably be identified 
with Spitzer features that were in later spectra. In particular, Spitzer emission lines features 
identified with corresponding ones in KIN inner spatial regime, dominated by light origi- 
nating in the close vicinity the WD, are well- represented by a cosmic distribution of atomic 
elements. We assumed that any condensates within the blast radius of the nova would be 
sublimated away and dissociated into atoms, and that any nucleosynthesis that occurred 
in the outer layers of the WD during TNR would negligibly impact abundances. The KIN 
inner spatial regime spectral features were keyed to Spitzer spectral features in the 4/16 and 
4/26 spectra. Similarly, the KIN outer spatial regime, light predomi nantly from a region ~ 



17 AU from the WD, is keyed to the Spitzer spectrum taken on 9/9 (lEvans et al.ll2007d ) as 
it is assumed that the abundances of atomic species in the latter, nebular Spitzer spectrum, 
would reasonably be representative of the environment observed by that KIN channel. 

Note that other symbiotic novae (e.g. V1016 Cyg, RR Tel) have broad emission features 
around 10 and 18 /im, evident in ou r Spitzer spectrum, t hat have bee n succ es sfully fitted 



by cry stalline silicate features (e.g. ISacuto et al.l (120071 ): ISchild et al.l (120011 ): lEyres et al. 



( 1l998n ). Motivated by this fact, we calculated the temperature and emission SED of various 
species of dust at the range from the pseudo-photosphere of the WD at which the KIN's 
outer spatial regime has maximum transmission. The center of the first constructive fringe, 
when the null fringe is located on top of the WD, at 9.8 fim and a projec ted baseline 



of 84 m, is ~ 12 mas from the WD. At a distance to the object of 1.4 kpc (IBarry et al. 



2008al ) this corresponds to about 17 AU. At day 3.8, the pseudo-photosphere of the WD 
has a luminosity of 1.6x10^ Lq (radius, 18.1 Rq, and temperature, 27050 K), computed 
by interpolating between values for the post-outburst evolution displayed in Table 5. The 
temperature of a black dust grain in thermal equilibrium is 272 K, using Eqn. A3 from the 
Appendices in this paper, and is well below the sublimation temperature of silicate dust, ~ 
1500 K. The dust feature in the Spitzer 9/9 spectrum is wide enough at FWHM that it falls 
across over seven spectroscopic elements in the KIN outer spatial channel meaning that the 
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dust is both spectrally resolved and spatially localized. 



Figure displays the association of identified Spitzer atomic emission lines with the 
continuum-normalized KIN inner and outer regime spectral features. Clearly the two traces 
are markedly different. Note that when the inner and outer spatial regimes are summed, the 
result closely follows the re-binned Spitzer spectra with the exception of the KIN outer spatial 
regime feature (lower trace) centered at 9.8 fim. Based on our assumptions of a primarily 
nebular environment in the vicinity of the outer spectrum, the atomic metals evident in the 
upper trace (KIN Inner spatial regime) would be unlikely to contribute to this feature. In 
any case, the total power in these metal lines in the wavelength range of this feature (cf. 
Table [3]), including Ca I, Si I, C I, & C II, evident in the upper trace is 0.05 Jy, while that 
in the broad feature centered at 9.8/im in the lower trace exhibits is 0.24 Jy. Our models 
suggest that the source of this feature may be hot silicate dust in the temperature range 800 
- 1000 K. 

Note that our KIN data detects the faint emission from well outside of the blast radius, 



assuming an initial shock front velocity of 3500 km/s (jO'Brien et al.l 120061 ) and negligible 
deceleration. The radiation from this spatial region originates primarily from material around 
the nova that has been illuminated and warmed by photons from the nova fiash and, as a 
result, must have existed before the nova event. This establishes that silicate dust, created in 
the vicinity of the RS Ophiuchi system some time previous to t he 2006 outburst, is detected 
by our measurements, and is consistent with the conclusions of lEvans et al.l (l2007d ). 



5. A physical model of the recurrent nova 

One aspect of the RS Oph binary system that has been neglected in the current literature 
regarding the recent outburst is the effect of the motion of the two stars through the wind 
created by mass loss from the red giant star. Garcia (1986) suggested that there could 
be a possible ring of material of diameter less than 40 AU around the RS Oph binary 
system or possibly surrounding the red giant component, based on his measurements of an 
absorption feature in the core of the Fe II emission line profile at 5197 A. The observations 
were performed in 1982 and 1983, several years before the 1985 nova outburst. 



Subsequently, and motivated by somewhat different observations, iMastrodemos fc Morris 



( 1l999l ) computed three-dimensional hydrodynamical models of morphologies of the envelopes 
of binaries with detached WD and RG/AGB components in general. Their purpose was to 
see if these models could reproduce some of the observed characteristics of axisymmetric or 
bipolar pre-planetary nebulae. Their study focused on a parameter space that encompassed 
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outflow velocities from 10 to 26 km/s, circular orbits with binary separations from 3.6 to 50 
AU, and binary companions having a mass range of 0.25 to 2 Mq. For binary separations 
of about 3.6 AU and mass ratios of 1.5, it was possible to generate a single spiral shock 
that winds 2-3 times around the binary before it dissipates at > 25 times the radius of the 
RG star. There is a density enhancement of about a factor of 100 over the normal density 
in the wind in the plane of the orbit of the two stars, and an under density or evacuated 
region perpendicul ar to the p l ane of the orbit. Observational support for this model was 



found recently by iBode et al.l (120071 ) who detected a double-ring stru cture in HST data. 



which they interpreted as due to an equatorial density enhancement and iMauron Sz Huggins 



(120061 ) who observed a spiral pattern around the AGB star AFGL 3068 both of which were 

del of iMastrodemos &: Mqrri^. T he underlying binary, a red giant 



consistent with the mo 
and white dwarf, was discovered by iMorris et al.l (120061 ) , who also determined the binary 



separation and hence approximate orbital period, whic h was consistent wit h expectations 



from the appearance of the spiral nebula pattern seen by lMauron fc Huggind and the model. 



Figure [6] shows the proposed geometry of the nebula in the plane of the o r bit of the 



RS Oph system based on the parameters adopted from iDobrzycka &: KenyonI (1l994j ) for 
the system, including an orbital period of 460 days, and an inclination angle of about 33 
degrees, for the epoch just prior to the nova outburst. The spiral shock model produces 
an archimedian spiral nebula, with the separation between adjacent windings of about 3.3 
mas based on the period noted above and a wind speed from the red giant star of 20 km/s. 
Based on these parameters we estimate that about 17 such rings could be created between 
outbursts, with the overall size of the nebula of the order of 100 mas. However, there could 
be as few as 10 rings to as many as about 20 rings, depending on the parameters, some of 
which are not well known. 

In order to better understand the Keck Nuller observations and other high angular 
resolution observations, we have modeled the evolution of the circumstellar nebula following 
the outburst, which we display in Fig. [91 For the purposes of this discussion, we adopt 
the model of Hachisu and Kato (2001), but with the simplification of a flat accretion disk 
geometry, i.e., not warped as in their paper. The details of our calculations are presented in 
the Appendices to this paper. 

The results are of fundamental importance to the interpretation of our KIN data. We 
begin with the most immediate effect after the blast, which is the sublimation of dust within 
a zone where the temperature of a blackbody dust grain would be > 1500 K in equilibrium. 
Figure [7] displays the evolution of the dust sublimation radius, which is roughly 4-5 mas 
until about day 70, after which it steadily declines to << 1 mas about 250 days after the 
outburst. This means that much or all of the dust within this zone has entered the gaseous 
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phase, except for a small "sliver" of material in the shadow of the red giant star. This 
provides additional hot gas (rich with metals) that is subsequently affected by the blast 
wave passing through within the next few days. 

As the two stars move relative to each other in their orbit about their common center 
of mass, the location of the shadow of the nova moves, and consequently the material in the 
shadow that has not been affected by the blast wave from the nova will be sublimated during 
this luminous phase, creating hot gas in the vicinity within a few mas of the stars. This 
material may still have the type of repeating density structure that was initially present, and 
may be observable with high angular resolution instrumentation. 

Whether or not this particular material is observable, depends in part on the density 
distribution as a function of latitude of the plane of the orbit of the two stars. If the 
material is uniformly distributed over Air steradians, then the fraction of the total solid 
angle subtended by the red giant star as seen by the white dwarf star is approximately 0.13- 
0.29% for a RG star between 27 and 40 R© and a binary separation of 1.72 AU. However, 
the hydrodynamical studies of Mastrodemos and Morris (1999) show that the density falls 
off steeply as a function of latitude, and is as low as 1% of the mid-plane density by latitudes 
of about ± 40 degrees. This gives a scale height of about 9 degrees. Thus, the red giant star 
subtends a much bigger fraction of the solid angle up to this scale height, i.e., to as much as 
1.9%. 

Furthermore, studies of the supernova blast waves around red giant stars indicate that 
the blast wave diffracts around the RG star, with a hole in the debris of angular size ~ 31-34 
degrees, which did not depend strongly on whether the companion star was indeed a RG 
star or a main sequence or subgiant star (Marietta, Burrows, & Fryxell 2000). In this case 
the fraction of the sky subtended by this hole in the debris field is ~ 10-13 %. If the material 
is concentrated in the mid-plane, the effect is substantially bigger as noted in the previous 
paragraph. Thus there are several reasons to expect that material in the shadow of the RG 
star can have an observable effect. 

Stripping of material from the red giant companion by the blast wave for supernova has 
been studied by Wheeler, Lecar, and McKee (1975) and by Lyne, Tuchman, and Wheeler 
(1992), however. Lane et al. (2007) showed that stripping of material for the RG star is 
negligible for the less energetic blasts from RS Oph. 

Another effect is the heating of the surface of the red giant star that faces the nova (see 
Appendices), as seen in Figure El for times past the maximum in the visible light curve. Our 
calculations indicate that this side of the red giant star increases from about 3400 K to as 
much as 4200 K within a few days past the outburst. The temperature steadily declines from 
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day 70 until it reaches equilibrium with the other side around day 250. The surface of the RG 
star is initially heated by the shock from the blast wave (not included in our calculations), 
however, the continued heating during the high luminosity phase is important as it affects the 
interpretation of data from modern complex instrumentation such as stellar interferometers 
and adaptive optics systems, where various subsystems are controlled at wavelengths other 
than the measurement wavelength. 

Figure M displays a schematic view of the system geometry from days 4 to 90 after the 
recent outburst. Assuming typical wind velocities of about 20 km/s for the red giant wind, 
there are roughly 17 rings separated by approximately 3.3 mas that form between RS Oph 
outbursts. The top left panel displays the system geometry at 4 days post-outburst. A gray 
ring is drawn in the center of the figure to indicate the size of the region affected by the 
blast at this epoch. In (a) the outer part of the spiral is overlaid with light gray to indicate 
that it is not known if the material stays in a coherent spiral past the first few turns. The 
diameter of the shocked region is about 8.8 mas assuming the blast wave tra vels at a velocity 



of < 1800 km/s in the plane of the orbit, using the velocity measured by IChesneau et al. 



(120071 ) ■ at approximately the same epoch as our measurements. (We acknowledge that most 
out-of-plane blast- wave velocities noted in the literature are much higher than this.) In this 
figure we assume the blast wave moves at constant velocity as it traverses the spiral shock 
material. We show the extreme case in which the blast wave is 100% efficient in sweeping 
up material, thus creating a ring-like structure that propagates outward from the system. 
Note this figure is meant only to be illustrative and differs in detai l from estimates of the 



position of the blast wave from observations, such as lO'Brien et al.l (120061 ) who obtained a 
value of the shock radius of 8.6 mas at day 13.8 by which time the blast wave had apparently 
slowed considerably. We obtain a value of 10.5 mas for the shock radius assuming a constant 
velocity of ~ 1800 km/s from one day past the initiation of the TNR process, which we tak e 



to be about 3 days before maximum light in V band (jStarrfield. Sparks fc TruranI Il985l ). 
The value of 8.6 mas at day 13.8 is consistent with a somewhat lower mean velocity, of the 
order of 1400 km/s. It is beyond the scope of this paper to compute the evolution of the 
blast wave, however, this figure makes a connection with the previous observational evidence 



for a "ring" of material concentrate d in the plane of t 



and with the recent observations of IChesneau et al.l (120071 ) who observe different velocities 



le or bit as discussed by Garcia (1986) 



perpendicular to the plane of the orbit than in the plane of the orbit. 

The evolution of the luminosity of the red giant and nova is also significant in aiding 
our understanding of the observations. Figure [TU] (a) displays computed V band light curves 
up to day 250 after the outburst including the effect of the irradiation of the accretion disk 
and the irradiation of the red giant star by the nova. Note this calculation overestimates the 
total luminosity of the nova during the period from about 10 days to about 50 days, and this 
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is likely due to the simplified disk geometry that we have employed in our own calculations. 
Most importantly these light curves show that the V band luminosity is dominated by the 
nova for about the first 50-70 days, and after that the red giant star dominates the V 
band luminosity. This is significant as most telescopes track on V band light (including 
interferometers) and the tracking center then moves by a mas or so during the post-outburst 
evolution of the system. H band luminosity evolution is plotted in Fig. [10] (b), and is 
different than that of the V band evolution. At H band the nova dominates the luminosity 
only for the first 5 days or so, and after that the H band light curve is completely dominated 
by the red giant star. This means the phase center for fringe detection is offset from that of 
the tracking center from days 5 onward by a mas or so, as mentioned above. The N band 
luminosity, displayed in Fig. [10] (c), evolves like that of the H band, and the nova dominates 
the mid-infrared light only up to day 4. After that there is an offset between the tracking 
center and N band fringe center like that noted for the H band. 



6. Implications to the Keck Nuller and other high-angular resolution 

observations 

We now reexamine data from the 2006 outburst within the framework of the spiral shock 
model of the geometry of RS Oph presented in the last section. In its simplest form, this 
geometry evolves into a bipolar morphology similar to that seen in the radio emission from 
the outburst that was discussed by O'Brien et al. (2006), as seen in Figure [9]^d). 

Such a geometry provides a natural explanation for some of the differences between the 
interferometric and other measurements, as the blast wave would be restricted and slowed in 
the orbital plane due to the high density regions, while its flow would be relatively unimpeded 
perpendicular to that plane. This corresponds well to what was observed by Chesneau et al. 
(2006), where the observed velocity components came from distinctly different position an- 
gles. These authors also calculate that the re d giant star should b e at position angle 150-170 



degrees at the time of the outburst. Recently, iBrandi et al.l (120081 ) have derived spectroscopic 



orbits for both components of the RS Oph system based on the radial velocities of the M 
giant absorption lines and the broad emission wings of Ha. They have also constrained the 
orbit inclination, i > 49 ± 3, using the estimated hot component mass, Mh sin^ i = 0.59 Mq, 
and assuming that the white dwarf mass cannot exceed the Chandrasekhar limit. Thus the 
plane of the orbit is such that the high velocity flow would be expected to be mostly East- 
West and the the continuum emission would have an elliptical shape with the position angle 
measured by Chesneau et al. (2006). Furthermore, the spiral shock wave tends to have the 
largest densities near the white dwarf star and on the opposite side of the red giant star. 
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the separation being of the order of two to four times the separation between the two stars, 
i.e., a few mas. The larger separation occurs when the red giant star has been spun up so 
that it is synchronous with the orbital motion. Thus it is possible that some of the data can 
be interpreted in terms of emission from hot clumps of material within the spiral shocks. 
Detailed radiative transfer calculations will need to be performed to refine this model. 

There is a wealth of recent observational work on RS Oph across the spectrum from X- 



rav (ISokoloski et al. 



2007 



20061 : iBode et al.ll2006 ) , near -infrared ( Evans et al.ll2007a : Chesneau et a 



Das et al.ll2006l ). radio ( O'Brien et al.ll2006l ). and now mid-infrared ( Evans et al 



2007b 



this work). Generally speaking, the observa tional picture is consistent with the shocked wind 
model as desc ribed bvlBode fc Kahru (119851) . Extensive modeling of the light curves has been 
performed by lHachisu &: Katd (l2000l . l200ll ). which include effects such as the irradiation of 
the red giant star by the white dwarf, and the accretion disk, which become important about 
4 days post-outburst, and which are the basis for our own calculations presented in Figs. 
I6|7t and [TOl In this general picture the recurrent nova differs from a classical nova in that 
the high-velocity ejecta from the WD is impeded by the wind from the red giant star, which 
in turn generates a shock wave that p ropagates t hrough the red giant wind. Observations of 
the shock wave in the near- infrared by lDas et al.l confirm this picture and trace the evolution 
of the widths of Pa P and O I lines as a function of time post-outburst, and clearly show 
the expected free-expansion phase of the shock ended about 4 days post-outburs t. This is 
also consistent with independent measurements conducted by iBode et al.l (120061 ) in which 
this phase was proposed to have ended by approximately day 6. 

Most of the observational work has used spectroscopic methods, but only the interfer- 
ometric observations are capable of spatially separating the various components of RS Oph 
that contribute to the emission that is seen spectroscopically. Chesneau et al. (2007) ob- 
served RS Oph 5.5 days post-outburst in the continuum at 2.13 /im, Br 7 at 2.17 /im, and 
He I at 2.06 /im with the AMBER instrument on the VLTI. They fitted their data with 
uniform ellipses, gaussian ellipses, and a uniform ring. The models had excellent consistency 
in terms of the position angle of the ellipse (~140 degrees) and the ratio between minor and 
major axes (0.6). For the continuum at 2.13 /^m, a uniform ellipse had a major axis of 4.9 
± 0.4 and minor axis 3.0 ± 0.3 mas, while the gaussian ellipse had a major axis of 3.1 ± 0.2 
and minor axis of 1.9 ± 0.3 (FWHM). These measurements are consistent with the expected 
size of the shocked region at the epoch of their measurements. 

Monnier et al. (2006) presented IOTA results in the near-infrared bands at H and K, 
and had somewhat different conclusions than some of the other workers because they were 
able to fit their visibility and closure phase data best with a binary model with two sources 
separated by 3.13 ± 0.12 mas, position angle of 36 ± 10 degrees, and brightness ratio of 
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0.42 ± 0.06. Their gaussian fits at 2.2 //m gave a FWHM of 2.56 ± 0.24 mas for the same 
period of time as the AMBER observations. A striking feature of those results is that the 
size of the emitting region decreased 10-20% between about days 4 and 65 post-outburst. 
For example, the size at 2.2 fim actually decreased from about 2.6 to 2.0 mas (FWHM), 
while at 1.65 fim, the size decreased from 3.3 to 2.9 mas. However the 2 fim continuum 
sizes are in approximate agreement. Monnier et al. (2006) rule out an expanding fireball 
model, however, they would have over-resolved the fireball anyway since it would be about 
8.8 mas in diameter, at a distance of 1600 pc, or a substantially larger angular size if the 
distance were smaller. Lane et al. (2007) observed RS Oph also with IOTA, PTI, and the 
Keck Interferometer at H band over a longer time period, up to 120 days from the V band 
maximum. They observed an increase in diameter from approximately 3 to 4 mas from days 
to 20 and a decrease to less than 2 mas in diameter around day 120. They interpreted the 
near-infrared size data in terms of a simplifed model of free-free emission in the postnova 
wind as the mass ejected in the wind decreased during that time period. 

The emitting region at 10 /xm is somewhat larger than that seen at 2 fim, for example, 
our data are fitted by a gaussian 4.0 ± 0.4 mas (FWHM), while the IOTA data had a size 
of 2.6 ± 0.2 mas. By 4 days post-outburst, a spherical shock wave would be expected to 
have a radius of about 3.9 mas assuming a speed of approximately 1730 km/s (O'Brien 
et al. 2006) and a distance of 1.4 kpc. Thus, according to this simple model, all the 
interferometric IR continuum emission should be coming from within the post-shock region. 
The AMBER/VLTI results indicate a more complex picture of the velocity field of the 
expanding material, with two indicated - a "slow" field between -1800 and 1800 km/s, and 
a "fast" one between -3000/-1800 km/s and 1800/3000 km/s. The position angle of the 
emitting material for the two velocity groups differ, with the "fast" component being well 
defined in the East-West (position angle 90 or 270 degrees ± 5 degrees) direction and the 
"slow" component with position angles from 55 to 110 degrees modulo 180 degrees. 

Recent observations by Bode et al. (2007) using the HST confirm the elliptical shape 
seen with AMBER, with an axial ratio of 0.5-0.6. The HST observations also confirm the 
velocity structure seen by AMBER, i.e., the "fast" velocity field in the East-West direction, 
of the order of 3000 km/s, and the "slow" velocity about half that of the "fast" one. The 
HST results showed that the expansion rate in the plane of orbit of the stars decreased, from 
about 0.62 mas/day on day 13.8 to 0.48 mas/day on day 155. This is a reduction in velocity 
from about 1700 km/s to about 1300 km/s. 

The spiral shock model is clearly relevant to the interpretation of the IOTA data pre- 
sented in Monnier et al. (2006). As noted above, the increased density of gas and dust in 
the arms of the spiral would certainly provide an impediment to the free expansion of the 
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fireball, as well as provide a reservoir of hot material that would emit strongly at 2 /im. 
Their data indicate a closure phase that is consistent with zero or 180 degrees for the first 
epoch (days 4-11) and convincingly non-zero only for hour angles from -1.5 to -1 hr on the 
second epoch (days 14-29), and at +1/2 hr on the third epoch (days 49-65). Hence, the 
closure phase signals and the binary interpretation may be more consistent with hot clumpy 
material that is cooling, some of which may (re-) condense into dust as the white dwarf's 
outburst luminosity declines from its most luminous state, i.e., as shown in Fig. [61 Note 
that the sublimation radius decreases from ~ 5 mas to about 2 mas from day 70 to day 
120. Furthermore, as the luminosity changes, the star tracker on IOTA may be providing a 
different optical center to the fringe detection system, since in the first few days during the 
outburst the optical emission is dominated by that centered on the white dwarf star itself, 
whereas after it has cooled, the optical emission is dominated by the red giant star. Thus 
the effects of changes in the optical tracking and interferometric phase center will need to 
be included in the analysis to properly interpret the data. The offset between the optical 
tracking center and the fringe phase center can cause a miscalibration of the visibility. Since 
the 2 /zm emitting region is actually decreasing in size with time, as seen in their Fig. 1 
and Table 2, it seems more likely they are observing the cooling of this hot material near 
the two stars than actually resolving the binary, however, the effects of the material in the 
shadow of the red giant star must be included in the interpretation of the near-infrared and 
mid-infrared data. 



7. Summary and Conclusions 

We analyzed data from the recurrent nova RS Oph for the epoch at ~ 4 days post- 
outburst using the new KIN instrument. These data allowed us to determine the size of the 
emitting region around the RS Oph at wavelengths from 8-12 fim. By fitting the unique 
KIN inner and outer spatial regime data, we obtained an angular size of the mid-infrared 
continuum of 6.2, 4.0, or 5.4 mas for a disk profile, gaussian profile (FWHM), and shell 
profile respectively. The data show evidence of enhanced neutral atomic hydrogen emission 
located in the inner spatial regime relative to the outer regime. There is also evidence of 
a 9.7 fj,m silicate feature seen outside of this region, which is consistent with dust that had 
condensed prior to the outburst, and which has not yet been disturbed by the blast wave 
from the nova. Our analysis of the observations, including the new ones presented in this 
paper, are most consistent with a new physical model of RS Oph, in which spiral shock waves 
associated with the motion of the two stars through the cool wind from the red giant create 
density enhancements within the plane of their orbital motion. 
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Purther observations are needed to clarify this new picture of the RS Oph system. One 
issue that has not been fully resolved is whether or not the red giant star really overflows 
its Roche lobe. If so a hot spot would be expected where the material streaming from the 
red giant envelope hits the accretion disk, and UV or X-ray observations could search for 
this effect. Another approach would be to observe RS Oph over several orbital cycles using 
infrared photometry to look for variations of the light curve due to the departure of the red 
giant star from spherical symmetry. High resolution spectra could also help. Confirmation 
of the rotational velocity of the red giant star measured by Zamanov et al. (2007) would be 
worthwhile, and could provide another estimate of the absolute size of the red giant star, 
assuming that it is co-rotating with the orbit. Further KIN observations within the next few 
years would also be helpful as they could show evidence of the re-establishment of the spiral 
shock wave, and perhaps some information about the shape of the circumstellar material and 
dust formation. Another epoch of HST observations would also determine the deceleration 
of the outflow in the two directions, i.e., within the plane of the orbit of the two stars and 
along the poles. 

Theoretical studies of the motion of the blast wave in an environment with a high density 
region in the plane of the orbit of the two stars are also worthwhile. In particular, it would 
be important to understand how the blast wave is diffracted around the red giant star and 
how it propagates in an medium with the periodic density enhancements in the plane due 
to the spiral pattern. 

The recurrent nova RS Ophiuchi is a rich system for the study of circumstellar matter 
under extreme physical conditions. Continued study will provide important insights into 
Type la supernovae, of which RS Oph may be a progenitor. 
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A. Luminosity Evolution 

In this appendix we discuss the mathematical formulation used to derive the evolution 
of the luminosity of the white dwarf star starting at the maximum magnitude in V band. 
Our calculations are based on those presented in the paper by Hachisu and Kato (2001), 
who expanded on their discussion of their light-curve model presented in Hachisu and Kato 
(2000). 

Our purpose is to evaluate the light curve not only at V band, but also at H, K, and 
N bands, which have been used by the AMBER, IOTA, Keck, and PTI interferometers to 
observe RS Oph after the 2006 outburst. Our discussion is in fact relevant to any modern 
instrument that has wavefront sensing and control, and also tracking, at wavelengths different 
than that of the observations of interest. For example, the IOTA interferometer performs 
its precision pointing at V band, but observations are made at H band. Similarly the Keck 
Interferometer does precision pointing at V band, wavefront control of the two large apertures 
at H band, and the data are taken at N band. This difference in wavelengths is important 
because as the V band luminosity of the white dwarf star decreases, the observed light from 
the binary is mainly from the red giant star, which gives a shift in the center of light. This 
shift can affect the calibration of these instruments as there is an offset, which will change 
over time as compared to a calibrator, which has all the same optical center at all times for 
all of these wavebands. 

The evolution of the stellar parameters for the white dwarf are given in Table [5] of 
Hachisu & Kato (2001). We assume the parameters for the red giant star remain constant 
with Rrg = AORq, and Trg = 3400 K. 

The blackbody flux density (erg cm~^ s~^ Hz~^) at a frequency, z/, for a star of temper- 
ature, T*, radius, R^, and distance, D is: 

F,(z/) = 271 {R J D)\huyc^)/[exp{hu/kT,) - 1], (Al) 

where h, k, and c, are Planck's constant, the Boltzmann constant, and the speed of light, 
respectively. 

The total luminosity (erg s~^) is: 

= 47rai?,'T/, (A2) 
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where a is Stefan's constant. 



A.l. Sublimation Radius 

The temperature of a black grain in radiative equilibrium at a distance, r, from a star 
with luminosity, L^, is: 

T \ T / R \ 



- ) - 71 V J ^^^^ 

A simple rearrangement of the above equation yields a formula for the sublimation radius, 
Rsub, assuming a sublimation temperature, Tsub- 

For this paper we assume a sublimation temperature of 1500 K. 

The evolution of the sublimation radius for the red giant star (dashed lines) and nova 
(solid lines) beginning at the maximum in V band is displayed in Fig. [71 Note that the 
sublimation radius remains approximately constant at about 5 mas for the first 70 days, 
after which it gradually reduces to < 0.2 mas after about day 120. 



A. 2. Illumination of the Red Giant Star 

During the high luminosity phase of the outburst, the surface of the red giant star facing 
the nova is heated substantially. We calculate the effect of the nova luminosity on the red 
giant star using Eqn. (10) of Hachisu & Kato (2001), which is based on thermal equilibrium 
between the faces of the two stars, which we display here: 

aTrgi"^ = Tjrg cos{9)L^d/ (47rr2) + aTrg'^ (A5) 

where T^gi is the new effective temperature of the red giant star for facing side and where 
we include an extra term from the white dwarf luminosity, L^^- There are two additional 
constants included in this new term, the first is rj^g, which is an effective emissivity for the 
stellar surface, while the second, cos(6') is the average inclination angle of the surface. The 
term, r, is the distance between the two stars, and T^g is the original temperature of the red 
giant star. Equation A5 can be rearranged into a simple form after substituting for L^^: 

T»/ + T,/ (A6) 



T 



rgt 



r]rg cos{9) 



R 



wd 
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Equation A5 was evaluated as a function of time, using the evolution of the white dwarf 
luminosity derived from Eq. A2 and the values from Table 5. The result of our computation 
of this effect is displayed in Fig. [HI in which we use r = 325/?©, cos(6') = 0.5, rjrg = 0.5, and 
Tr„ = 3400 K. 



A. 3. Accretion Luminosity 



The computation of the evolution of the accretion luminosity is based on the treatment 
of Hachisu & Kato (2001) with the simplification of a flat accretion disk, instead of adding 
the complexity of a warped accretion disk that is used i n the ir paper. Thes e calculations 
are based on the well known iLynden-Bell fc Pringld (119741 ) and lPringle I (jl98ll ) equations for 
the luminosity of an accretion disk, where the accretion luminosity is based on a numerical 
integration of the flux from the accretion disk, assuming a temperature distribution across 
the disk. Normally the temperature of the accretion disk is determined solely by the accretion 
rate onto the disk for normal viscous heating of the disk. To this Hachisu & Kato (2001) 
added a term based on additional heating of the disk from the white dwarf star as it is in 
a high luminosity state and is contracting and heating during the constant high luminosity 
period after the outburst. The flux density of the disk at frequency u is given by: 



Fdisk 



27r cosf 



Rou 



pBy[Tdisk{p)]dp 



(A7) 



The conventional expression for the radial distribution of the temperature of the disk is: 



Tdiskl (p) 



/ 3GM^M 



P 



1/2 



1/4 



(A^ 



where Tdiskl is the temperature in the accretion disk (jPringle Ill98ll ) and D is the distance 
to the Earth. The quantity B^, is the Planck function, M* is the mass of the white dwarf, R^, 
is its radius, M is the accretion rate onto the disk, and a and G are the Stefan-Boltzmann 
constant and the gravitational constant, respectively. The maximum temperature in the 
accretion disk occurs at a radius Rmax = 1-36 R^,, and is given by 



T 

J- rr 



0.488 



3GM,M 



1/4 



(A9) 



as discussed in iLynden-Bell &: Pringld ( 119741 ) , iPringle I (Il98ll ) , and iHartmann I (119981 ) . In 
this treatment we neglect the emission from the boundary layer, which occurs over a very 
small angular region around the surface of the white dwarf star and emits at a very high 
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temperature, which would have a neghgible contribution to the luminosity in the infrared 
region near 10 /xm. The evolution of the luminosity of the accretion disk will be computed 
using the formulation of Hachisu and Kato (2001), as set up below. 

Let us now formulate the evolution of the accretion disk luminosity by adjusting some 
of the parameters in the equation above as a function of time. In particular the outer radius 
of the accretion disk Rdisk is parametrized by a power law decrease from day 7 until day 79 
using the following formulae: 

Rdisk = aRi* (AlO) 



a 



f 



(t-to)/72 



(All) 



In these equations R\* is the inner critical radius of the Roche lobe for the white dwarf 
(nova) component, which we assume is 138.6 Rq The parameter a helps define the size of 
the accretion disk which varies between = 0.1 and ai = 0.008, and so the power law form 
of Eqn. AlO is to make an interpolation function. Haschisu and Kato include irradiation of 
the accretion disk as the WD luminosity changes between days 7 and 79, where the accretion 
disk temperature is given by: 



0-[Tdisk{p)t = Vdisk COs{9inc) 



L 



wd 



3 GM,M 
87rp3 



R* 
P 



1/2 



(A12) 



The second term of this equation is the traditional accretion luminosity term, while the first 
term includes the luminosity of the white dwarf. The parameter rjdisk is an efficiency factor, 
and is assumed to be 0.5, and cos(^inc) is 0.1, based on the average inclination of the surface. 
They assume the outer radius of the accretion disk is at a temperature of 2000 K, and is 
not affected by radiation from the WD photosphere. For the remainder of the calculation 
we will compute a in discrete time intervals corresponding to what we have done before over 
the days that this effect matters. 

The integral is numerically evaluated at the center wavelengths for V (0.55 /im), H (1.65 
fim), and N (10.5 /im) bands, using an inner radius given by the radius of the white dwarf 
star and the outer radius given by the equation for Rdisk, for the time steps of Table 5 and 
the white dwarf parameters in that table. The results arc plotted in Fig. 10. In this figure 
the quantities plotted include the the contributions of the red giant and white dwarf alone 
(stars and sohd squares), the red giant including effects of irradiation by the white dwarf 
(diamonds), the accretion disk alone (sohd circles and dash-dot-dot line), the irradiated 
accretion disk (triangles and dashed line), and the total for the white dwarf including the 
irradiated accretion disk (diamonds and solid line). 
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Table 1. Observing Log for RS Ophiuchi. 



Object 


Type 


Time 


U 


V 


Airmass 






(UT) 


(m) 


(m) 




Chi UMa 


cal 


15:07:07 


23.85 


80.61 


1.38 


Chi UMa 


cal 


15:15:39 


21.92 


81.25 


1.41 


RS Oph 


trg 


15:50:15 


54.57 


64.75 


1.46 


RS Oph 


trg 


16:03:46 


55.35 


64.37 


1.39 


RS Oph 


trg 


16: 12:35 


55.75 


64.12 


1.35 


Rho Boo 


cal 


16:34: 24 


39.63 


75.14 


1.08 



Table 2. RS Ophiuchi model fitting results. 



Source Model Angular Size (mas) Radiant Flux Major Size (mas) Minor size*^ (mas) 

(Nband) (Jy) (K band) (K band) 



Uniform Disk 6.2 ± 0.6 22.4 ± 3.9 4.9 ± 0.4 3.0 ± 0.3 

Uniform Gaussian^ 4.0 ± 0.4 22.4 ± 3.8 3.1 ± 0.2 1.9 ± 0.3 

Uniform ShelF 5.4 ± 0.6 22.4 ± 3.8 3.7 ± 0.3 1.9 ± 0.2 



'^Sizes for continuum values at 2.3 /im after IChesneau et al 
^FuU width at half maximum. 

■^Spherical shell with thickness 1.0 mas - optically thin. 



(12007t ). 
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Table 3. Mid Infrared Spitzer Line List: N-band. 



Wavelength 


ID 


Spitzer 4/16 

(Jy) 


Spitzer 4/26 

(Jy) 


Spitzer 9/9 

(Jy) 


7 /ifin 


Pfn— HT- (\ — ^ 
1 J LX II 1 . yj o 


LiLyljLyL.ljLyU. 


H o"t" or*"!" oH 

U.L.ljLyL.ljLyLi 


LLLy LLyL. LLyLl 


7 fi52 


\NpVT] 

i V C- r -i 


r]p-(-pp-(- pri 


np-j-ppf pn 


Lie tCL^ tCLl 


s isn 

O. -LOU 


FpT 


n nnq 


n nn7 




8 7fin 

o. / uu 


ffT- 10 — 7 

11 1 . ±u / 


n 049 


u.uoo 


090 ^ 


8 Q8.'i 


\VTT] Si T }ir Ca Tl 






010 


Q ni 7 




n ni ^ 


n 091 




Q 

y . z,oo 


Cn T 


n nn4 


n nn4 

u.uu^ 




Q 407 


O LI 


n ni4 


n ni 

U .U -LO 








n nriQ 


n nns 

u.uuo 




Q 790 


(ITT 

1 1 


n m 8 


n ni 7 










00^ 




10.285 


Mo/ 


0.007 


0.006 




1 n 4Q9 


ff r- 1 9 _ c 

1 . ±Z( o 






01 




NpT 

1 V C--£ 


096 


096 




10.833 


C7] 


0.018 


0.016 





11.284 


if/: 9 - 7 






0.007 


11.318 


Hel 


0.060 


0.053 




11.535 


Hell 


0.016 


0.013 




12.168 


HI:36 - 11 


0.011 


0.009 




12.372 


Hua -HI:7-6 


0.158 


0.150 


0.037 


12.557 


Sil 


0.040 


0.036 




12.803 


[Nell] 






0.159 


12.824 


Hel 


0.009 


0.015 




13.128 


Hell 


0.009 


0.011 




13.188 


HI: 18 - 10 


0.009 


0.010 


0.005 



'^Detected in Spitzer spectrum but not fit due to intrusion of band edge. 
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''Blended with neutral hydrogen hnes at 8.721 and 8.665 iim 

Note. — Some of these species were blended with others, and, due to the 
difficulty in deblending Spitzer low-resolution channel spectra, may in some cases 
be misidentified. This is not critical, however, to the conclusions of this paper. 
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Table 4. Continuum-normalized, Mid-infrared KIN Emission Source Identification. 



Center Wavelength Spectroscopic Width Flux Attributed to 

(/xm) (/xm, FWZM) (Jy) 



KIN Inner S]:)a.tial Regimo 



8.7 


8.3 - 9.1^ 


0.06*^ 


H I: 10-7, Fel, Ca P 


9.4 


8.9 - 11.1 


0.02 


Ca I, Si I, C I, C II 


10.4 


9.9 - 11.1 


0.04 


Mg I, Ne I, C I] 


11.4 


11.1 - 11.8 


0.07 


He I, He II 


12.5 


12.1 - 12.6 


0.15 


Hua 


KIN Outer Spatial Regime 


8.9 


8.3 - 9.5 


0.14 


HI: 10-7, [VII], Si I, Ca I]*^ 


9.8 


9.0 - 10.7 


0.24 


Silicate Dust'' 


11.4 


11.0 - 11.8 


0.19 


H I: 9 - 7, He I, He II 



'^Approximate FWZI continuum crossing points. 

^Because Keck and Spitzer data were not taken simultaneously and because 
of the extreme transient nature of the RN outburst it is not possible to assert 
particular flux numbers to the atomic lines attributed to the KIN feature. The 
flux listed should be considered the maximum, bounding flux for any one of the 
atomic species shown here. 

■^All atomic line emission in the KIN inner spatial regime is assumed to be 
emitted by a species at approximate cosmic abundance as for RN only a very 
moderate amount of nucleosynthesis is theorized to occur. 

"^AU atomic line emission in the KIN Outer spatial regime are assumed to be 
predominantly of nebular abundance with some contribution from uncondensed 
metals. 



-31 - 



•^Spectral feature is spectrally and at least partially spatially resolved. Silicate 
dust feature has more than seven spectral elements across it and emits relatively 
strongly at a distance centered ~ 17 AU from the WD. All flux is assumed to be 
attributable to emission by the dust. 
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Table 5. Post-outburst Evolution of the White Dwarf Star in RS Oph 



Timc'^ 


Radius 


Temperature 


(days) 


(Ro) 


(K) 





45 


12200 


2 


21 


20000 


14 


1.6 


67000 


29 


0.56 


114000 


51 


0.23 


181000 


72 


0.083 


302000 


119 


0.0037 


1020000 


251 


0.003 


350000 



^Stellar parameters for the white dwarf star after the 
2006 outburst of RS Oph based on Hachisu & Kato 
2001 

'^The zero time for these parameters is the V band 
maximum hght, which occurs about 3-4 days after the 
thermonuclear runaway process begins. 
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Keck Interferometer Nuller 




Fig. 1. — The upper part of plot displays a conceptual view of the operation of the KIN 
system. There are two short baselines, separated by 5 m, and two long baselines separated by 
85 m. A vr phase shift is applied to the pupils on the long baseline. The autocorrelation of the 
four pupils is shown as well as the equivalent intensity pattern, also called the transmission 
pattern. Radiation passing through the white striped regions is detected by the system. 
Lower part of the plot displays the sequence used for the measurement process. There are 
two chopping sequences, (a) and (b) display the ~ 5 Hz chopping between the null and 
bright fringe patterns on the source using the 85 m baseline for starlight subtraction, used 
to measure the null response, while (c) and (d) show the ~1 Hz chop sequence between 
the two short basehnes used to remove the telescope and sky backgrounds. Note that the 
long-baseline fringes are not to scale with the field of view in this depiction. 
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Fig. 2. — Plot of three best- fit continuum models against Keck NuUer data. All three models; 
disk, shell, and Gaussian, that minimized simultaneously against the null leakage and 
intensity spectra, lie effectively on top of one another. The top trace gives the dimensionless 
null leakage or interferometric observable which is the null fringe output divided by the 
intensity spectrum. The lower trace is the constructive fringe output or intensity spectrum. 
As described in the text we have removed several data points associated with emission 
features from the original set for the purpose of fitting the continuum. 
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Nova RS Ophiuchi Mid Infrared Spectra 




Fig. 3. — Plot of KIN spectra together with a Spitzer spectTum. Upper trace is the intensity 
spectrum or inner KIN spatial regime - light dominated by the inner 25 mas about the 
center of the source brightness distribution at mid band. The middle trace is the nulled 
fringe output which is the interferometric observable times the intensity spectrum multiplied 
by 2 (to correct for transmission through the fringe pattern for extended sources) to give 
the source brightness in the outer spatial regime. This is predominantly emission from 
material greater than about 12. 5 mas (^17 AU at 1.4 kpc) from the center of the source. 
The lower trace is Spitzer data (lEvans et al.ll2007bl ) from day 63 boxcar averaged to yield 
approximately the same spectral resolution as the Keck Interferometer NuUer. None of the 
data were continuum normalized. Note the features between 9 and 11 fim and how the 
inner and outer spatial regime spectra are different from one another. Note that the Spitzer 
spectra have very low spatial resolution and combine the light from the entire region detected 
by both KIN inner and outer spatial regimes. 
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X 



Spitzer Observations 
RS Ophiuchi 





J I I I I \ I I I I 

10 12 14 

Wavelength (nnicrons) 



Sept 9 



1 6 



Fig. 4. — Spitzer space telescope spectra from days 63, 73 and 209. Here we see that the 
continuum drops rapidly as time advances with the spectral emission features almost identical 
on days 4/16 and 4/26. Data obtained on 9/9 is starkly di f ferent with a strong solid-state 
feature evident between 9 and 13 microns (see lEvans et al. (|2007ch for more details on the 
Spitzer data). 
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Nova RS Ophiuchi KIN Mid Infrared Spectra 




Fig. 5. — Continuum normalized Keck Interferometer Nuller spectra 3.8 days after peak 
V-band brightness. Spectra are shown matched with hues identified in Spitzer spectra. The 
inner spatial regime line identification was matched to the cosmic abundances seen in the 
earlier Spitzer spectra while the outer spatial regime was biased towards the later spectrum 
on September 9, 2006 that included a solid state feature and were assumed to have primarily 
nebular abundances. Please see the text for more discussion of the fine identification process. 
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RS Ophiuchi 
Before Nova 




H 100 mas »l 

Fig. 6. — Proposed model of circumstellar material surrounding the binary RS Oph before 
the nova eruption. The interaction of the white dwarf and red giant star in the slow dense 
wind of the red giant star can create a spiral shock with an enhanced density in the plane 
of the orbit of the two stars. The overall size of the dense in-plane material is of the order 
of 100 mas if RS Oph is at a distance of 1400 pc and the wind speed is about 20 km/s. 
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Fig. 7. — Sublimation radius of dust (mas) as a function of the number of days post-outburst 
for the Nova and for the red giant companion. 




Fig. 8. — Temperature of red giant star for the surface that is facing the Nova as a function 
of the number of days post-outburst. Also plotted is the temperature of the non-irradiated 
side of the red giant. Modeled temperature rise is due to heating from irradiation by the 
nova and neglects any effect due to the passage of the forward shock. 
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RS Ophiuchi 
Day 4 




■ 100 mas 



RS Ophiuchi 
Day 57 




RS Ophiuchi 
Day 21 




' 100 mas 



RS Ophiuchi 
Day 90 




100 mas 



4 



■ 100 mas ■ 



Fig. 9. — Spiral shock wave model of RS Oph. (a) Top left panel displays the system 
geometry at 4 days post-outburst. A gray ring is drawn in the center of the figure to 
indicate the size of the shocked region at this epoch. The outer part of the spiral is overlayed 
with light gray to indicate that it is not known if the material stays in a coherent spiral past 
the first several turns. The diameter of the shocked region is about 5 mas assuming the blast 
wave travels at a velocity of 1730 km/s in the plane of the orbit, (b) Top right panel. The 
blast wave is now about 13 mas in diameter on day 21. (c) Bottom left panel. The blast 
wave is about 36 mas in diameter on day 57. (d) Bottom right panel. By day 90 the blast 
wave has traversed the entire spiral pattern. 
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Fig. 10. — (a) V Band luminosity as a function of time for Nova and red giant components 
of binary system. (b) H band luminosity evolution, (c) N band luminosity evolution. 



